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We study the possibility to stabilize unstable steady states and unstable periodic orbits 
in chaotic fractional-order dynamical systems by the time-delayed feedback method. By 
performing a linear stability analysis, we establish the parameter ranges for successful sta- 
^1 . bilization of unstable equilibria in the plane parametrizad by the feedback gain and the 

time delay. An insight into the control mechanism is gained by analyzing the characteristic 
^ . equation of the controlled system, showing that the control scheme fails to control unsta- 

ble equilibria having an odd number of positive real eigenvalues. We demonstrate that the 
^«0 , method can also stabilize unstable periodic orbits for a suitable choice of the feedback gain, 

providing that the time delay is chosen to coincide with the period of the target orbit. In 
r^ , addition, it is shown numerically that delayed feedback control with a sinusoidally modu- 

Qh' lated time delay significantly enlarges the stability region of the steady states in comparison 

^ . to the classical time-delayed feedback scheme with a constant delay. 

bJQl PACS numbers: 05.45. Gg, 02.30. Ks, 45.10.Hj 

o 

The fractional order dynamical systems have attracted remarkable attention in the last decade. 
CN ' Many authors have studied the chaotic and hyperchaotic dynamics of various fractional order 

systems, such as those of Duffing, Lorenz, Rossler, Chua, Lii, Chen, etc., which are introduced 
Q^ \ by changing the time derivative in the corresponding ODE systems, usually with the fractional 

OO ' derivative in the Caputo or Riemann-Liouville sense of order < a < 1 [l|-l5!]- One interesting 

problem is to analyze the lowest value of parameter a under which fractional order dynamical 
systems show chaotic or hyperchaotic behaviors. Stability analysis, synchronization and control of 
fractional order systems by using different techniques are also widely investigated [QtIJI and are of 
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INTRODUCTION 



great interest due to their application in control theory, signal processing, complex networks, etc 

30]. 

KA \ In this paper we investigate the possibility to control unstable equilibria and unstable periodic 

^^ ■ orbits in fractional order chaotic systems by a time-delayed feedback method. Pyragas introduced 

^ the time-delayed feedback control (TDFC) in 1992 by constructing a control force in a form of 

a continuous feedback proportional to the difference between the present and an earlier value 
of an appropriate system variable r/, i.e. K[r]{t — T) — 7?(t)], with K and T being the constant 



control parameters denoting the feedback gain and the time delay, respectively |21l423l|. To stabilize 
unstable periodic orbits, the delay T is chosen to match the period of the unstable orbit. In the 
case of controlling unstable equilibria, the optimal delay is related to the intrinsic characteristic 
time scale given by the eigenfrequencies of the uncontrolled system. In both cases, the control 
force vanishes when the target state is reached, rendering the method noninvasive. In this sense. 
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the unstable states of the uncontrolled system are not changed, since the control force acts only if 
the system deviates from the state to be stabilized. 

The main advantage of the Pyragas method over the other control methods is that it does not 
require the knowledge of the system's equations and the positions of the unstable states whose 
control is desireable. The method has been used in many concrete applications 24l433l|. and also 
rigorously investigated analytically 3J-l50f| . In recent works 5l|, l52| , it has been shown that the 
efficiency of the method is greatly improved by deterministic or stochastic modulation of the time 
delay T. This variable delay feedback control (VDFC) has been shown successful in stabilization 
of unstable equilibria in systems described by both ordinary and delay differential equations. 

This paper is organized as follows. In Sec. II we present the time-delayed feedback controller for 
stabilization of unstable equilibria in fractional-order dynamical systems. We first look at a general 
non-diagonal case of such a control scheme and derive the stability conditions for the equilibrium 
points in the presence and absence of control. To illustrate the control method, we choose a 
fractional order Rossler system in a chatoic regime controlled via a single state variable with a 
Pyragas-type feedback force. By performing a linear stability analysis of the controlled system, 
we calculate the domain of successful control of the unstable equilibria in the plane parametrized 
by the feedback gain and the delay time. It is shown, both numerically and analytically, that the 
unstable equilibrium point with an odd number of positive real eigenvalues cannot be controlled 
by the time-delayed feedback method for any values of the control parameters, thus extending 
the validity of the odd-number limitation theorem 4a, |47t] to the case of fractional-order systems. 
In Section III we perform a numerical calculation of the control domain by the Pyragas delayed 
feedback with a modulated time delay in a form of a sine-wave. In Section IV we give a numerical 
evidence that a time-delayed feedback control can successfully stabilize unstable periodic orbits, 
and estimate the corresponding feedback gain intervals for which such a control is possible. A 
summary of the obtained results is given in Section V. 



II. DELAYED FEEDBACK CONTROL OF UNSTABLE EQUILIBRIA 

A. Stability analysis 

We consider a general n-dimensional nonlinear fractional-order dynamical system under a non- 
diagonal form of delayed feedback control in the sense of Pyragas: 



f D^^xiit) =/i(x(t))+Fi(t), 
^D^-Xn{t) =/„(x(i)) + F„(t), 



(1) 



where 



Fiit) = Y,Kij[xj{t-T)-x,{t)] 



(2) 



is the delayed feedback force applied to the ith component of the system, consisting of contributions 
of all the system components, Kij are the gain factors of the feedback terms, T is the constant time 
delay, x = (xi, X2, . . . , x„) is the state vector, and f = (/i, /2, . . . , fn) is the nonlinear vector field 
that determines the dynamics of the unperturbed system. The notation D" is the time fractional 



derivative in the Caputo sense defined as [5. 
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where m = [a], i. e., m is the first integer which is not less than a. In the following, we consider 
the fractional orders Oj to be in the interval (0, 1). 

In the case Kij = K6ij, where 5ij is the Kronecker delta, the generalized control scheme ([1]) 
reduces to TDFC with a diagonal coupling, and when the control force is applied only to a single 
system component and consists only of contributions of the same component, it yields the original 
TDFC control scheme introduced by Pyragas. 

Let P = {xl,X2, . . . ,x*) be an arbitrary equilibrium point of the system ([1]) in the absence of 
control (Kij = 0), being a solution to the nonlinear algebraic system: 



' /l(xi,X2,... ,Xn) = 0, 
f2{xi,X2,...,Xn) =0, 

_ fn{xi,X2,... ,Xn) =0. 



(6) 



Assuming that P is an unstable equilibrium point of the uncontrolled system, we wish to find 
the domain in the parameter space of the feedback gains Kij and the time-delay T for which P 
becomes locally asymptotically stable under TDFC force ([2]). The stability of P under a non- 
diagonal feedback control ([IJ-(I2]) can be determined by linearizing ([T]) around P, which leads to 
the linear autonomous system: 
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is the Jacobian matrix of the free-running system, with aij = {dfi/dxj) calculated at P, Xi(t) = 
Xi{t) — X* are the transformed coordinates in which the equilibrium point is at the origin, and Fi{t) 



are the components of the feedback control force in the new coordinates, i.e. 

n 

m) = J2Kij[x,{t-T)-x,{t)]. 



(9) 



By applying the Laplace transform to Eqs. ([7]) and by using the formula for the Laplace transform 
of the fractional derivative in the Caputo sense [5Q |: 

m—l 

C[D2xi{t)] = s'^X.is) - J2 xf (0+)s"-l-^ (10) 

A;=0 

where Xi{s) = C[xi(t)] are the Laplace images and x\ (0) are the initial conditions, we obtain: 

A(s)-X(s) = B(s), (11) 

where 



A(.) 



/s"i-aii+Kii(l-e-^^) -ai2 + Ki2(l-e-^^) ... -ai„ + Ki„(l - e"^^) \ 
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(12) 
represents a characteristic matrix of system ([7]), X(s) = col(Xi(s), ^2(5), . . . , X„(s)) is the column 
vector of the Laplace images, and 



B(.) 



X2(0)s"2-i + e-^^ Y.%1 K2j JIt x,{t)e~''dt 



(13) 



For the sake of the argument, let us assume that all the roots s of det[A(s)] = are positioned in 
the left complex s-plane (Re(s) < 0). The assumption is equivalent to the claim that det[A(s)] / 
for Re(s) > 0, meaning that A(s) is invertible matrix in the right complex s-plane, and thus 
X(s) = A(s)^^ • B(s) has a unique solution X(s) for Re(s) > 0. Furthermore, it is easy to check 
that 



from which we have 



lim s'Xis) = lim [A(s)"^ ■ sB(s)l = 0, 



lim x(t) = lim sX(s) = 0. 



(14) 



(15) 



by using the final-value theorem for the Laplace transform. Thus, we prove the following theorem. 
Theorem 1. The equilibrium point P of the system ([I])-([2]) is locally asymptotically stable if 
and only if all the roots s of the characteristic equation: 



det [A(s)] = 0, 



(16) 



have negative real parts, i.e. 

|arg(s)|>7r/2. (17) 

The matrix A(s) is given by Eq. (fT2]) . and its components aij = (dfi/dxj) are evaluated at the 
equilibrium point P. D 

Remark 1. For the case when there is no control {Kij = 0), the conditions for stability of a 
particular equilibrium point is still provided by Theorem 1. In this case, the characteristic matrix 
of the uncontrolled system simplifies to 

A(s) = S-I-A, (18) 

where S = col(s"^, s°^, . . . , s""), I is n x n identity matrix, and A is the Jacobian matrix of the 
unperturbed system evaluated at the equilibrium point, and given by Eq. dS]). In the special case 
cti = Oi2 = ■ ■ ■ = Oin = Oi, the characteristic Eq. (jlSp can be recast into the form 



det f AI - Aj = 0, (19) 

which is a polynomial of degree n in A^where A = s". In this case, the stability condition (J17p can 
be rewritten in the form of Matignon ja| 

|arg(A)| >a7r/2, (20) 

meaning that the stability region is bounded by a cone, with vertex at the origin, extending into 
the right half of the complex A-plane such that it encloses an angle of iba7r/2 with the positive 
real axis. Therefore, the equilibrium point of the unperturbed system is stable if and only if all 
the roots of the characteristic polynomial (J19p are placed outside this cone. 

To derive another important result related to the limitation of the time-delayed feedback method 
dH)-©, we consider the function det[A(s)], with A(s) given by Eq. ([T2|), and s G M+. We take P 
to be an unstable equilibrium of ([T]) in the absence of external perturbation {Kij = 0), and A the 
Jacobian matrix of P. One can easily deduce that 

lim det[A(s)] = +oo, (21) 

s— >+oo 

and 

n 

lim det[A(s)] = det[-A] = TT(-ei), (22) 

s-i-0+ ■*■ , 

1=1 

where e^ are the eigenvalues of A. Evidently, if A has an odd number of positive real eigenvalues, 
then lim3_^o+ det[A(s)] < 0. In this case, the sign of det[A(s)] is changed from negative to positive 
when s sweeps the real interval [0,-|-cxd). Since det[A(s)] is a smooth function in s, there exists 
at least one positive real root of the characteristic equation det[A(s)] = 0, meaning that the 
equilibrium point P cannot be stabilized by the time-delayed feedback controller ([2]). The result 
is summarized in the following theorem. 

Theorem 2. (odd-number limitation) Let P be an unstable equilibrium point of the 
fractional-order system ([T]) in the absence of control {Kij = 0) , and A the corresponding Jacobian 
matrix at P. If A has an odd number of positive real eigenvalues, then the time-delayed feedback 
control ([2]) cannot stabilize the unstable equilibrium P for any values of the control parameters 
Kij and T. D 



The result is an extension of the odd-number limitation theorem 46|, |47| to fractional-order 



systems with respect to unstable fixed points. We note that the odd-number limitation has recently 



been refuted by Fiedler et al. 48|, |49[] for the case of unstable periodic orbits in systems described 



by ordinary differental equations. 
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FIG. 1: The phase plots of the chaotic attractor in the free-running fractional-order Rossler system. The 
parameters are a = 0.4, b = 0.2, c = 10 and a — 0.9. 

B. Numerical example 

To illustrate the time-delayed feedback control in fractional order chaotic systems, we consider 
a fractional order Rossler system in the form: 



D^xit) =-y{t)-z{t), 

D^y{t) =x{t) + ay{t)+F{t), 
D^zit) = z{t) [x{t) - c] + 6, 



(23) 



where 



F{t)=K[y{t-T)-y{t)] 



(24) 



is the Pyragas feedback controller applied through a single component (y-channel), and a, b and c 
are the parameters of the free-running system. In the following, we take a = 0.4, b = 0.2, c = 10 
and a = 0.9, for which the uncontrolled system {K = 0) has a chaotic attractor (see Fig. I) [3|]. 

The unperturbed Rossler system has two equilibrium points Pi = (x^,y!j_,z^) and P2 = 
{x*_,y*_^z*_), where 
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Linearization around the equilibrium points leads to the following linear autonomous system: 




(28) 




where 



(29) 

is the Jacobian matrix, and x{t) = x{t) —Xj_, y[t) = y{t) —y±, z{t) = z{t) — Zj_ are the transformed 
coordinates in which the corresponding fixed point is at the origin. According to Eq. ()20p . the 
equihbrium point of the hnearized system (j28p is asymptoticahy stable if and only if | arg(A)| > 
aTr/2 for all the eigenvalues A of the Jacobian matrix A. The equilibrium point Pi = (x!j_, y!j_, z!^) = 
(9.9919,-24.98,24.98) has eigenvalues Ai = 0.3844 and A2,3 = 0.0038 ± 5.0964i. It is an unstable 
saddle point of index 1 since Ai > 0, |arg(A2,3)| = 1.5701 > air/l. The equilibrium point P2 = 
{x*_,y*_,z*_) = (0.008,-0.02,0.02) has eigenvalues Ai = -9.9900 and A2,3 = 0.1990 ± 0.9797*, and 



it is an unstable saddle point of index 2 since Ai < 0, | arg(A2,3)| = 1.3704 < aTr/2 [551]. 

In the presence of TDFC, the linearized version of the system (I23p around the equilibrium points 
states: 

D^y{t) = A ■ y{t) \+K\ y{t - T) - y{t) | . (30) 

D-zit)J \z{t)J \ 



According to Theorem 1, the zero solution of system (I30p is asymptotically stable if and only if all 
the roots s of the characteristic equation: 

det [A(s)] = (31) 

have negative real parts, i.e. | arg(s)| > vr/2, where the characteristic matrix A(s) is given by: 
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A(s) = I -1 s'^-a + i^(l-e-*^) |. (32) 

c) 



The characteristic Eq. (I3ip can be numerically analyzed to obtain the domains of control for the 
unstable steady states Pi ^2 in the plane parametrized by the feedback gain K and the time delay 
T. 

In the absence of control, the equilibrium point Pi has an odd number (one) of positive real 
eigenvalues, and according to Theorem 2, it cannot be stabilized by the TDFC method. This 
result has been confirmed by a numerical analysis of the characteristic Eq. ()3ip , showing absence 
of stability domain in the {K, T) parameter plane. This observation is further confirmed by a 
numerical simulation of the system (I23p under TDFC (I24p . On the other hand, the fixed point P2 
can be controlled by TDFC, and the resulting stability domain is shown in Fig. 2. The stability 
islands (shaded areas) denote the values of the control parameters K and T for which all the 
eigenvalues s of the characteristic Eq. (j3ip are lying on the left complex s-plane, thus satisfying 
the stability condition (I17p . For these values of the control parameters, the control of the fixed point 
P2 is successful. As a verification, we performed a computer simulation of TDFC by numerically 
integrating the system (123p - (|24p . The resulting diagrams are shown in Fig. 3. The simulations 
were done by using a predictor-corrector Adams-Bashford-Moulton numerical scheme for solving 
fractional order differential equations [SSj]. Panels (a), (b) and (c) depict the dynamics of the state 
variables x(t), y{t) and z{t), respectively, and panel (d) shows the corresponding time series of the 
control signal F{t). In the simulations, the control parameters were K = 2 and T = 3, belonging 
to the domain of successful TDFC control depicted in Fig. 2. As expected, the simulation confirms 




FIG. 2: The stability domain in the {K^ T) parametric plane of the unstable equilibrium P^ in the fractional- 
order Rossler system (^51) under time-delayed feedback control ([M)) . The shaded areas correspond to the 
control parameter values for which stabilization of the fixed point P2 is achievable. The parameters of the 
free-running system are a = 0.4, h = 0.2, c = 10 and a = 0.9. 



a successful stabilization of the unstable equilibrium P2- Moreover, as indicated from panel (d) 
in Fig. 3, the control signal F{t) vanishes when the control is achieved, meaning that the control 
scheme is noninvasive. 

We note that the above analysis has been repeated for different parameter values of the free- 
running system. In each case, the resulting stability domains computed from Eqs. (|3ip -()32 p are in 
agreement with the numerical simulation of the TDFC method. Specifically, for a = 0.4, h = 0.2, 
c = 10 and variable a, we observed a decrease in the stability region as a is increased from a = 0.9 
to a = 1. On the other hand, as a becomes smaller than 0.9, the complex-conjugate eigenvalues 
of the equilibrium point P2 eventually escape the instability region described by the Matignon 
formula ()20p . resulting in a stable equilibrium P2 even without control. The critical value a = ac 
that corresponds to this eigenvalue-crossing of the conic surface between the different stability 
regions can be calculated from Eq. (|20p . In this case, Oc = 0.8724. 



III. VARIABLE-DELAY FEEDBACK CONTROL OF UNSTABLE EQUILIBRIA 



In recent papers [51 



52l ]. it has been demonstrated, both numerically and analytically, that the 
original Pyragas TDFC scheme can be improved significantly by modulating the time delay in an e 
interval around some nominal delay value Tq. In both deterministic and stochastic variants of such 
a delay variation, the stability domain was considerably changed, resulting in an extension of the 
stability area in the control parameter space if appropriate modulation is chosen. In the following, 
we will demonstrate numerically the successfulness of this variable-delay feedback control in the 
case of fractional-order chaotic Rossler system (I23p with a = 0.4, h = 0.2, c = 10 and a = 0.9. In 
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FIG. 3: Successful stabilization of the unstable equilibrium P2 in the fractional-order Rossler system by a 
time-delayed feedback control, (a)-(c) Time plots of the state variables x{t), y{t) and z{t). (d) The feedback 
control force F{t) vanishes when the control is achieved. The uncontrolled fractional-order Rossler system 
is in a chaotic regime, with parameters: a — 0.4, b — 0.2, c — 10 and a — 0.9. The control parameters were: 
K = 2, T — 3. The control was activated at t = 50. The total time span shown in each panel is 250 time 
units. 



this case, the feedback force F{t) is given by: 

F{t)=K[y{t-Tit))-y{t)], 
where we choose a time-varyng delay T(t) in a form 

T{t) = To + esm{ujt), 



(33) 



(34) 



modulated around a nominal delay value Tq with a sine-wave modulation of amplitude e and 
frequency oj. Obviously, if e = then T{t) = Tq = const, and the variable-delay feedback control 
is reduced to the classical Pyragas TDFC scheme. With this choice of the feedback force F(t), the 
control parameters of the proposed variable-delay scheme are K, Tq, e and w, and thus, the control 
parameter space is four-dimensional. For visualisation purposes, we may fix two of the control 
parameters and investigate the stability domains in the parametric plane spanned by the remaining 
two control parameters. To demonstrate the superiority of variable-delay feedback control over 
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FIG. 4: Domain of successful variable-delay feedback control in the {K, Tq) plane for the unstable equilibrium 
P2 in the chaotic fractional-order Rossler system (a = 0.4, b — 0.2, c = 10, a = 0.9). The control delay T{t) 
is modulated with a sine-wave, with e = 1 and w = 10. Combinations of K and Tq where the control scheme 
successfully stabilizes the fixed point P2 are plotted in grey. Note the shift of the origin along the Tq axis 
by an amount equal to e due to the physical limitation Tq > e of the controller. 



TDFC, we fix the modulation amplitude £ and the frequency uj and investigate the control domain in 
{K, Tq) parameter space. Numerical simulations show that the stability area is gradually increasing 
as e is increased from zero. Figure 4 shows such a control domain for e = 1 and a; = 10. The 
grey region indicates those values of the control parameters K and Tq for which the control of the 
unstable equilibrium P2 is successful. The stability domain is obtained by numerically integrating 
the linearized system 



D^xit) 
D^t) 




+ K{ y{t 





T{t)) 




m 



(35) 



with the Jacobian matrix A given by Eq. (|29p . It is evident that the control domain is significantly 
enlarged in comparison to the one in TDFC in Fig. 2. We note that numerical integration of the 
system ()35p for the unstable equilibrium Pi shows failure of the variable-delay control scheme for 
any values of K and Tq, suggesting validity of the odd- number limitation theorem also in the case 
of a time- varying delay. 

As a demonstration of the variable-delay feedback control in the fractional-order Rossler system, 
in panels (a)-(c) of Fig. 5 we show the dynamics of the state variables x{t), y{t) and z{t) for K = 3 
and Tq = 7, fixing the modulation amplitude e = 2 and frequency a; = 10. The time series indicate a 
successful stabilization of the unstable equilibrium P2- The method of control is again noninvasive, 
as indicated by the vanishing feedback force F{t) in panel (d) of Fig. 5. We note that for these 
parameter values, the control via TDFC (e = 0) is unsuccessful, as can be perceived from the 
stability domain in the TDFC case depicted in Fig. 2. 
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FIG. 5: A simulation of the variable-delay feedback control in the chaotic fractional-order Rossler system 
(a = 0.4, b — 0.2, c = 10, a = 0.9). The delay modulation is in a form of a sine- wave with amplitude e = 2 
and frequency uj = 10. The time series of the variables a;(i), y{t) and z{t) depicted in panels (a)-(c) indicate 
a successful control of the unstable equilibrium P2- The vanishing feedback force F{t) in panel (d) shows 
the noninvasiveness of the control method. The control parameters were: K = i, T^ ~ 7 . The control was 
activated at i = 70. The total time span shown in each panel is 200 time units. 



IV. DELAYED FEEDBACK CONTROL OF UNSTABLE PERIODIC ORBITS 



The Pyragas delayed feedback control method was originally aimed to stabilize unstable periodic 
orbits embedded into the chaotic attractor of the free-running system 2^. For this purpose, the 
time delay in the feedback loop was chosen to coincide with the period of the target orbit. By 
tuning the feedback gain to an appropriate value, the stabilization is achieved and the controller 
perturbation vanishes, leaving the target orbit and its period unaltered. 

In this section, we will give a brief demonstration of the Pyragas method to control unstable 
periodic orbits in the fractional-order Rossler system (J23p -( f24l) . As in the previous discussion, we 
use a = 0.4, h = 0.2, c = 10 and a = 0.9, for which the system is chaotic in the absence of external 
perturbation. 

In order to estimate the periods of the unstable orbits which are typically not known a priori, we 
use the fact that the signal difference F{t) = y{t — t) — y{t) at a successful control asymptotically 
tends to zero if the delay r of the controller is adjusted to match the period T of the target orbit. 
The method consists of calculating the dispersion (-7^^) of the control signal at a fixed value of 
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FIG. 6: Example of the "spectroscopic" procedure to determine the periods of unstable periodic orbits and 
their control domains under TDFC in the fractional-order Rossler system (a = 0.4, b = 0.2, c = 10, a = 0.9). 
(a) The dependence of the dispersion of the control signal T{t) = y{t — t) — y{t) upon the time delay t 
for K = 0.2 reveals segments of finite length (FP) related to a fixed point stabilization (compare with the 
T-intervals in Fig. 2 at i^ = 0.2), and a resonance peaks at those values of t corresponding to the periods 
of unstable periodic orbits. The period-three resonance point becomes more pronounced for smaller values 
of K. (b) Calculation of the stability domain for the period-one orbit. Note the peak a,t K ^ 0.195, for 
which the control of the orbit is the most robust. 



the feedback gain K for a given range of values of the delay r, excluding the transient period 
22, ISTJ. The resulting logarithmic plots of the dependence of the dispersion (^^/ on the delay r 
may contain several segments of finite r-intervals for which (^^) is practically zero, and a sequence 
of isolated resonance peaks with very deep minima. The former correspond to the stability domain 
of the fixed point P2, and the latter are the points at which r coincides with some accuracy to the 
periods of the unstable periodic orbits in the original system. The estimated values of the periods 
T can be made more accurate if one repeats this "spectroscopy" procedure for a larger sampling 
resolution of the r interval encompassing the resonance peaks. In this way, we have obtained the 
periods of the unstable period-one, period-two, and period-three orbits: Ti « 6.2, T2 ~ 12.49 and 
T3 ~ 18.89. The plot of the dispersion (J~'^) vs the delay r for K = 0.2 is shown in panel (a) of 
Fig. 6. 

The same approach could be used to calculate the intervals of the feedback gain K for which 
the corresponding orbits can be stabilized with the Pyragas controller. In panel (b) of Fig. 6 
we depict the dependence of the dispersion {J^"^) on the feedback gain K when the delay time r 
coincides with the period of the first unstable periodic orbit r = Ti = 6.2. In this case, the interval 
of the parameter K for which the orbit is stabilized is estimated to he K = [0.14,0.65]. A similar 
analysis yields K = [0.1,0.22] for a period-two, and K = [0.06,0.15] for a period-three orbit. It is 
observed that the control interval of the feedback gain K becomes narrower as the period of the 
target orbit is increased. 
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In Fig. 7 we show the results of the stabilization of period-one orbit (T = 6.2) for K = 0.25. 
Panels (a) and (b) show the projection of the system trajectory in xy and xz planes, respectively, 
after the control of the target period-one orbit has been established. The time-series of the state 
variables are given in panels (c)"(e), and panel (f) shows the feedback force that vanish after the 
controller is switched-on, warranting a noninvasiveness of the control procedure. Analogous results 
related to stabilization of period-two and period-three orbits are given in Figs. 8 and 9. 

A detailed bifurcation analysis of the chaotic Rossler system described by ordinary differential 



equations and subjected to a time-delayed feedback control has been performed recently 571], 
revealing multistability and a large variety of different attractors that are not present in the free- 
running system. A similar analysis in the case of fractional-order chaotic systems is left for future 
studies. 

V. SUMMARY AND CONCLUSIONS 

We have shown that the time-delayed feedback control can be used to stabilize unstable steady 
states and unstable periodic orbits in fractional-order chaotic systems. Although the control 
method was illustrated specifically for the fractional order Rossler system, it has also success- 
fully been applied to stabilize unstable equilibria and unstable periodic orbits in various other 
fractional-order dynamical systems. In all the cases, delayed feedback control with a variable time- 
delay significantly enlarges the stability region of the steady states in comparison to the classical 
Pyragas TDFC scheme with a constant delay. 

We find that equilibrium points that have an odd number of positive real eigenvalues cannot 
be stabilized by TDFC for any values of the feedback control parameters. The result is known 
as the odd-number limitation theorem, which extends to the case of fractional-order systems, as 
purported by Theorem 2. The odd-number limitation is also confirmed numerically in the case of 
a variable-delay feedback control. 

An analytical treatment of delayed feedback control of unstable periodic orbits in fractional- 
order systems is still lacking, and constitutes a promising subject for a future research. Applying 



the extended versions of the delayed feedback controller [37|, l39| to fractional-order systems is 
another interesting topic not tackled in this paper. This is especially important regarding the 
observations for the system used in this paper, that the control domains are becoming smaller 
for higher orbits, such that the periodic orbits of periods higher than three practically cannot be 
stabilized by the original controller. 

A detailed analysis of the variable-delay feedback control in fractional-order systems, including 
a theoretical understanding of the method and numerical computation of the stability domains 
in different parameter planes and for different types of delay modulations are also left for future 
studies. 
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FIG. 7: A successful stabilization of period-one orbit in the chaotic fractional-order Rossler system by a 
time-delayed feedback control. The control parameters are K = 0.25 and T = 6.2. Panels (a)-(b) depict 
the phase plots of the controlled periodic orbits in the post-transient regime. Panels (c)-(e) depict the time 
series of the state variables x{t), y{t) and z{t). The vanishing feedback force F{t) in panel (f) demonstrates 
the noninvasiveness of the control method. The parameters of the unperturbed system are a = 0.4, h ~ 0.2, 
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time units. 



15 



[9] W. Deng, C. Li and J. Lu, Nonlinear Dyn. 48, 409 (2007). 

[10] W. Deng and C. Li, J. Phys. Soc. Jpn. 74, 1645 (2005). 

[11] W. Deng and C. Li, Physica A 353, 61 (2005). 

[12] W. Deng, Phys. Rev. E 75, 056201 (2007). 

[13] T. Zhou and C. Li, Physica D 212, 111 (2005). 

[14] C. Li, W. Deng and D. Xu, Physica A 360, 171 (2006). 

[15] R.S. Barbosa, J.A.T Machado and I.M. Ferreira, Nonhnear Dyn. 38, 305 (2004). 

[16] R.S. Barbosa, J.A.T Machado and A.M. Galhano, J. Vibration and Control 13, 1407 (2007). 

[17] E.J.S. Pires, J.A.T Machado and P.B. Oliveira, Signal Processing 83, 2377 (2003). 

[18] C. Li, W. Sun and J. Kurths, Physica A 361, 24 (2006). 

[19] X. Wang, Y. He and M. Wang, Nonhnear Analysis 71, 6126 (2009). 

[20] Y. Tang, Z. Wang and J. Fang, Chaos 19, 013112 (2009). 

[21] E. SchoU and H. G. Schuster (ed.). Handbook of chaos control (Wiley- VCH, Weinheim, 2008), second 
completely revised and enlarged edition. 

[22] K. Pyragas, Phys. Lett. A 170, 421 (1992). 

[23] K. Pyragas, Phil. Trans. R. Soc. A 364, 2309 (2006). 

[24] K. Pyragas and A. Tamasevicius, Phys. Lett. A 180, 99 (1993). 

[25] S. Bielawski, D. Derozier, and P. Glorieux, Phys. Rev. E 49, R971 (1994). 

[26] T. Pierre, G. Bonhomme, and A. Atipo, Phys. Rev. Lett. 76, 2290 (1996). 

[27] K. Hall, D. J. Christini, M. Tremblay, J. J. CoUins, L. Glass, and J. Billette, Phys. Rev. Lett. 78, 4518 
(1997). 

[28] D.W. Sukow, M. E. Bleich, D. J. Gauthier, and J. E. S. Socolar, Chaos 7, 560 (1997). 

[29] J. M. Krodkiewski and J. S. Faragher, J. Sound Vib. 234, 591 (2000). 

[30] T. Fukuyama, H. Shirahama, and Y. Kawai, Phys. Plasmas 9, 4525 (2002). 

[31] M. G. Rosenblum and A. S. Pikovsky, Phys. Rev. Lett. 92, 114102 (2004). 

[32] C. von Loewenich, H. Benner, and W. Just, Phys. Rev. Lett. 93, 174101 (2004). 

[33] O. V. Popovych, C. Hauptmann, and P. A. Tass, Phys. Rev. Lett. 94, 164102 (2005). 

[34] P. Hovel and E. Scholl, Phys. Rev. E 72, 046203 (2005). 

[35] S. Yanchuk, M. Wolfrum, P. Hovel and E. SchoU, Phys. Rev. E 74, 026201 (2006). 

[36] T. Dahms, P. Hovel and E. Scholl, Phys. Rev. E 76, 056201 (2007). 

[37] J. E. S. Socolar, D. W. Sukow and D. J. Gauthier, Phys. Rev E 50, 3245 (1994) 

[38] K. Pyragas, Phys. Lett. A 206, 323 (1995). 

[39] A. Ahlborn and U. Parlitz, Phys. Rev. Lett. 93, 264101 (2004). 

[40] A. Ahlborn and U. Parlitz, Phys. Rev. E 72, 016206 (2005). 

[41] H. G. Schuster and M. B. Stemmler, Phys. Rev. E 56, 6410 (1997). 

[42] K. Pyragas, Phys. Rev. Lett. 86, 2265 (2001). 

[43] K. Pyragas, V. Pyragas, I. Z. Kiss and J. L. Hudson, Phys. Rev. Lett. 89, 244103 (2002). 

[44] K. Pyragas, V. Pyragas, I. Z. Kiss and J. L. Hudson, Phys. Rev. E 70, 026215 (2004). 

[45] W. Just, T. Bernard, M. Ostheimer, E. Reibold and H. Benner, Phys. Rev. Lett. 78, 203 (1997). 

[46] H. Nakajima, Phys. Lett. A 232, 207 (1997). 

[47] H. Nakajima and Y. Ueda, Physica D 111, 143 (1998). 

[48] B. Fiedler, V. Flunkert, M. Georgi, P. Hovel and E. Scholl, Phys. Rev. Lett 98, 114101 (2007). 

[49] W. Just, B. Fiedler, M. Georgi, V. Flunkert, P. Hovel and E. Scholl, Phys. Rev. E 76, 026210 (2007). 

[50] C. M. Postlethwaite and M. Silber, Phys. Rev. E 76, 056214 (2007). 

[51] A. Gjurchinovski and V. Urumov, Europhys. Lett. 84, 40013 (2008). 

[52] A. Gjurchinovski and V. Urumov, Phys. Rev. E 81, 016209 (2010). 

[53] M. Gaputo, Elasticita Dissipacione (Zanichelli, Bologna, 1969). 

[54] K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equa- 
tions (Wiley, New York, 1993). 

[55] W. Deng and J. Lli, Chaos 16, 043120 (2006). 

[56] I. Podlubny, Fractional Differential Equations (Academic Press, San Diego, 1999). 

[57] A. G. Balanov, N. B. Janson and E. SchoU, Phys. Rev. E 71, 016222 (2005). 

[58] K. Diethelm, N. J. Ford and A. D. Freed, Nonlinear Dyn. 29, 3 (2002); K. Diethelm, N. J. Ford, A. D. 
Freed and Y. Luchko, Comput. Methods Appl. Mech. Engrg. 194, 743 (2005). 



16 




-15-10-5 5 10 15 

x(t) 
15 



15-10-5 5 10 15 

JC(t) 




FIG. 8: Results of stabilization of period-two orbit in the chaotic fractional-order Rossler system. The 
control parameters are K = 0.12 and T = 12.49. The parameters of the unperturbed system are as in Fig. 
7. The control was activated at t = 50. The total time span shown in each panel is 250 time units. Note 
that one of the period-two peaks in panel (e) is barely visible. 
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FIG. 9: Results of stabilization of poriod-thrcc orbit in the chaotic fractional-order Rossler system. The 
control parameters are K = 0.08 and T = 18.89. The parameters of the unperturbed system are as in Fig. 
7. The control was activated at t = 50. The total time span shown in each panel is 300 time units. Note 
that one of the period-three peaks in panel (e) is too small to be seen on the scale of this figure. 



